Pair condensation and bound states in fermionic systems 
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We study the finite temperature-density phase diagram of an attractive fermionic system that 
supports two-body (dimer) and three-body (trimer) bound states in free space. Using interactions 
characteristic for nuclear systems, we obtain the critical temperature Tc2 for the superfluid phase 
transition and the limiting temperature T^s for the extinction of trimers. The phase diagram features 
a Cooper-pair condensate in the high-density, low-temperature domain which, with decreasing den- 
sity, crosses over to a Bose condensate of strongly bound dimers. The high-temperature, low-density 
domain is populated by trimers whose binding energy decreases toward the density-temperature do- 
main occupied by the superfluid and vanishes at a critical temperature Tcs > Tc2- 



I. INTRODUCTION 

Pairing correlations and three-body bound states are 
universal properties of attractive fermions that are of con- 
siderable interest in a number of fields. Recent progress 
achieved in trapping and manipulating cold fermionic 
atoms has opened a new window on the many-body prop- 
erties of dilute Fermi systems 0. The possibility of ma- 
nipulating the strength of the interactions in these system 
by tuning a Feshbach resonance allows one to explore the 
phase diagram in different regimes, and in particular the 
crossover from strong to weak coupling in a controlled 
manner 

The properties of dilute nuclear matter play key roles 
in supernova and neutron-star physics, especially through 
its composition and equation of state at relevant temper- 
atures and pressures An understanding of few-body 
correlations in infinite nuclear matter is also important 
to the study of less tractable finite systems, notably nu- 
clei far from the valley of beta stability . In this work 
we shall focus attention on phenomena arising from two 
types of correlations in dilute isospin-symmetric nuclear 
matter, namely the formation of a condensate of Cooper 
pairs due to "pairing correlations," and the appearance 
of bound three-body clusters. Our primary goal is to 
identify the regions in the temperature-density phase di- 
agram where pair condensation, two-body bound states 
(dimers), and three-body bound states (trimers) are im- 
portant. The actual phase diagram of dilute nuclear mat- 
ter is likely to be much more complex, since in addi- 
tion one must anticipate the formation, in the nuclear 
medium, of bound alpha-particle clusters (mass number 
yl = 4) and higher- A clusters 0. However, the physics of 
pairing and three-body clusters and their interplay is of 
interest in its own right, so we shall not consider cluster- 
ing with A > i. Another problem domain in which pair- 
ing and bound-state formation are of interest is the high- 
density QCD of deconfined quark matter, which could be 
realized in the interiors of massive compact stars @ . 

Since the binding force is different in the diverse sys- 
tems mentioned above, the dependence of the binding en- 
ergy on the number of bound fermions is non-universal. 
For example, in nuclear systems the binding energy in- 



creases with the number of nucleons A to A — A. This 
is followed by a gap in the binding energy at ^ = 5, fol- 
lowed by another gap at ^ = 8. In contrast, for the QCD 
problem the three-quark states dominate the low-energy 
limit of the theory, while higher quark-number states are 
extremely rare. 

In the cold, dilute systems of fermionic atoms, a three- 
body bound state can be created if three different atoms 
or three different hyperfine states of the same atom are 
trapped. 

Partial- wave analysis of the nucleon-nucleon scattering 
data yields information on the dominant pairing channel 
in nuclear- and neutron-matter problems in a given range 
of density. At high densities, corresponding to labora- 
tory energies above 250 MeV, the most attractive pairing 
channel is the tensor-coupled 3P2 — 3F2 channel J^^ ^ , 
provided the isospin-symmetry is slightly broken; for per- 
fectly symmetric systems the most attractive pairing in- 
teraction is in the 3-D2-wave 0|. At low density, isospin- 
symmetric nuclear matter exhibits pairing due to the at- 
tractive interaction in the 3S'1-3I?1 partial wave [lol |. a 
tensor component of the force again being responsible for 
the coupling of S and D waves. This interaction channel 
is distinguished by the fact that it supports a two-body 
bound state in free space - the deuteron. Since a Cooper 
pair at finite chemical potential carries the same quantum 
numbers as the deuteron in free space, one anticipates 
a crossover from BCS pairing of neutrons and protons 
at high densities to clustering into deuterons and their 
Bose-Einstein condensation at low densities [Tol UlL IT^ . 
The neutron-proton condensate is fragile in dense nuclear 
matter, since the ratio of the gap to the chemical poten- 
tial is small (of order 0.1), and a small isospin asymme- 
try, refiected in a mismatch of the chemical potentials of 
neutrons and protons, disrupts the coherent superfiuid 
state 0, While pairing still exists in other chan- 
nels (15*0 at low and 3P2 at high densities), there are 
no bound states associated with these channels in free 
space 15]. The evidence for neutron-proton pairing in 
finite nuclei is seen in the excess binding of nuclei with 
N ^ Z 16]. 

Going one step further in the hierarchy of clusters re- 
quires treating three-body bound states in the nuclear 



2 



medium. As is well known, the non-relativistic three- 
body problem admits exact solution in free space [r^ll8l| . 
Faddeev equations sum the perturbation series to all or- 
ders with a driving term corresponding to the two-body 
scattering T-matrix embedded in the Hilbert space of 
three-body states. The counterparts of these equations 
in many-body theory were first formulated by Bethe [T9l | 
to gain access to the three-hole-line contributions to the 
nucleon self-energy and the binding of nuclear matter. 
In this approach, the Brueckner G-matrix is employed as 
the driving term in the three-body equations |2d| . More 
recently, alternative forms of the three-body equations 
in a background medium have been developed that ei- 
ther (i) use an alternative driving for ce ( the particle-hole 
interaction or scattering T-matrix) 123. l2^ or/and 
(ii) adopt an alternative version of the free-space three- 
body equations, known as the Alt-Grassberger-Sandhas 
form j2J, I23 • Our initial task will be to derive the homo- 
geneous integral equations that determine the in-matter 
bound-state wave-function and the corresponding eigen- 
states using the real-time Green's functions formulation 
of the in-matter three-body equations |23 . 

The paper is organized as follows. In Sec. ^ we 
study pairing in the isospin-singlet spin-triplet state in 
dilute finite-temperature nuclear matter. In Sec. Illll the 
three-body Faddeev-type equations for the bound-state 
problem at finite temperature and density are obtained. 
These equations are solved for a Malfliet-Tjon [2^ po- 
tential, and the density and temperature dependence of 



the bound-state energy and the three-body wave func- 
tion are explored. Section ll V I combines and summarizes 
the results of Sees, m and IIIII in a phase diagram of di- 
lute nuclear matter that supports pair correlations and 
three-body bound states. 



II. PAIRING 

This section begins with a brief description of the the- 
ory of nuclear superfluidity at finite temperature. Our 
aim here is to clarify the approximations entering the 
equations to be solved numerically. Readers familiar with 
the formalism can proceed to Subsec. Ill Bl for specifica- 
tion of the interactions adopted in computations, and to 
Subsec. FlI CI for the numerical methods, the results, and 
their analysis. 



A. Formalism 

We shall work within the real-time Green's function 
formalism, in which the propagators are assumed to be 
ordered on the Schwinger-Keldysh real-time contour [23| ■ 
Such an ordering is equivalent to arranging the correla- 
tion functions and self-energies in 2 x 2 matrices. The 
one-body Green's function matrix is defined in terms of 
the fermionic fields ijj{x) as 



c.„,Mx, x)-y ^>^^^^ x')J-[ {^l{x)^l{x')) {T-^Ux)M^')) J ' 

where T'^/°- are the time-ordering and anti-ordering operators, (...) stands for statistical averaging over the equilibrium 
grand-canonical ensemble and x is the space-time four-vector, while the Greek indices stand for discrete quantum 
numbers (spin, isospin). In equilibrium, the physical properties of the system are described by the retarded propagator 

G^^{x,x')^e{t-t')[G>^{x,x')-G<f,{x,x')] , (2) 

where 9{t) is the step function. (The retarded and advanced propagators are also related to the elements of the 
Schwinger-Keldysh matrix through a rotation in the matrix space by the unitary operator U = {1 + iay)/^/2, 
where ay is the y component of the vector of Pauli matrices.) The one-body propagator in the superfluid state is a 
2x2 matrix in Gor'kov space, 

f Gc.0ix,x') F^p{x,x')\ _( -i{Ti^^{x)i,l{x')) {iP^{x)i>p{x')) \ 
U[x,x) y_Fl^(^^^x') G^,{x,x')J [ mx)4{x')) -^{fMM(^')) ) ^ 

where Gapix, x') and F'^pix, x') are referred to as the normal and anomalous propagators. The 4x4 matrix Green's 
function satisfies the familiar Dyson equation 

G„^(x,x') = G°^(x,x') + ^ jd'x''d'x'''G^^{x,x''')^.^,{x''' ,x'')G,p{x'' ,x'), (4) 

7,(5 

where the free propagators G°^(a;, x') are diagonal in the Gor'kov space; the underline indicates that the propagators 
and self-energies are matrices in this space. We are restricting considerations to uniform fermionic systems, so that the 
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propagators depend only on the difference of their arguments by translational symmetry. A Fourier transformation of 
Eq. Q with respect to the difference of the space arguments of the two-point correlation functions leads to on- and 
off-diagonal Dyson equations 

Gafiip) = Goaisip) + Goajip)['Sjsip)Gsf3ip) + A^sip)Fdsi3ip)] , (5) 
FKf3ip) = Goa.ji-p)[A\sip)Gsf3{p) + 'Sys{-p)FAsfj{p)], (6) 

where p is the four-momentum, Goa/3(p) is the free normal propagator, and Sq,^(p) and Aq^(p) are the normal and 
anomalous self-energies. Summation over repeated indices is understood. The Dyson equations for the components 
Gapip) and Fajsip) follow from Eqs. © and ® through the time-reversal operation. Specifying the self-energies in 
terms of the propagators closes the set of equations consisting of (0 and © and their time-reversed counterparts. 
Before doing so, we can find the quasiparticle excitation spectrum of the superconducting phase in terms of yet 
unspecified self-energies. The spectrum is determined by the poles of the retarded propagators G^p{p) and F^p{p). 
Below we consider spin triplet, isospin singlet (neutron-proton) pairing in the tensor ZSl — ZDl channel, which implies 
Ft = 1 ® iTy F^ , where 1 is a unit matrix in spin space and iTy is the component of the Pauli matrix in the isospin 
space. For a spin-isospin conserving interaction, G^p = SapG and = (^a/jS. The solutions of Eqs. (jSJ and © in 
the quasiparticle approximation, which keeps only the pole part of the propagators, are 

= ul {uj - uj± + if^)-^ + vl{uj ~ + ii^)-^ , (7) 
F^ = F^ = UpVp [{lj -LJ+ + irj)-^ -{LJ-LJ-+ irj)-^] , (8) 
Gf = vl{uj-uj±+ iri)-^ + ul{uj -uj^+ (9) 



when expressed in terms of the quasiparticle spectrum lli± = ±y/E{p)'^ + A'^{p) and the Bogolyubov amplitudes Up 
and Vp normalized by = 1/2 + E{p)/2uj+ and Up + Vp — 1. The advanced propagators follow from the retarded ones 
via the replacement -\-iri — > —irj. In equilibrium, the elements of the 2x2 matrix appearing in Eq. ^ are determined 
from the retarded and advanced propagators as 

G<{p) = [G^ip) - G«(p)]/(a,), G>{p) = G«(p) + G<(p), (10) 
G^(p) = G<(p) + G^(p), G^{p)^G<{p)~G^{p), (11) 

where /(w) is the Fermi distribution function. For time-local interactions, both the pairing interaction (which we 
shall approximate by a two-body potential V{p,p')) and the pairing gap are energy independent. 

I 



The mean-field approximation to the anomalous self- 
energy (the gap function) is then 

A'^ip) - 2 J ^V{p,p')lmF^{cj,p')f{u;). (12) 

Substituting Eq. (jHl, we arrive at the quasiparticle gap 
equation. Further progress requires partial-wave expan- 
sion of the interaction. We keep the interaction in the 
coupled 35'1 — 3D1 channels to obtain two coupled inte- 
gral equations for the gap {I — 0, 2), 

MP) - - j %i:y.f^\p^p') 

(13) 

where D'^{k) = (3/87r)[Ag(fc) + A|(fc)] is the angle- 
averaged neutron-proton gap function and V^^^^{p^p') 
is the interaction in the 35'1 — 3Z?1 channel (the domi- 
nant attractive channel in dilute and isospin-symmetric 



nuclear matter). Below we shall work at constant tem- 
perature and density. The chemical potential is then de- 
termined self-consistently from the gap equation (|13|l and 
the expression for the density, 

n = -8/ ^ImG?(c.,p)/(.) 

= V ^[ulf{u:+)+vlf{u^.)\. (14) 

The factor 4 comes from the sum over the two projections 
of spin and of isospin. 

B. Interactions 

The gap equation and the three-body bound states 
have been studied for the Malfliet-Tjon (MF) potentials, 
whose simple form facilitates numerical solution of the 
three-body equations. These potentials fit basic proper- 
ties of few-nucleon systems, including the S'-wave phase 
shifts and the binding energies of the deuteron and triton. 
The MF potentials, central and local in r-space, consist 
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of a sum of attractive and repulsive Yukawa potentials, 



i=l 



Their p-space Fourier transform for 
has the form 



(15) 



partial waves 



Vip,p') 



TTpp' 



i=l 



Specifically, the driving term V^^^lp^p') in the gap equa- 
tion was taken as the MF-III parameterization of the 
potential itT^ . which reproduces the phase-shifts in the 
SSI channel and the binding energy of the deuteron, 
Ed = -2.23 MeV. (Since the MF potentials do not con- 
tain a tensor component, 5 and D states are not coupled, 
and the deuteron quadrupole moment is not reproduced.) 
The MF-III parameter values are: gi = —3.22, Ai = 1.55 
fm-i, g2 = 7.39, A2 = 3.11 frn'^ We also used the 
MF-V parametrization, which differs from MF-III in the 
strength of its attractive interaction, gi — —2.93. 



C. Solving the gap equation 

A number of algorithms exist for solution of the gap 
equation. If the potential is in separable form (or, being 
local, is approximated by a separable form), the origi- 
nal integral equation reduces to a set of algebraic equa- 
tions (^]. Alternatives available for nonseparable poten- 
tials include the eigenvalue method @|, the separation 
method '2^ , and the method of successive iterations [33| ■ 
Direct iteration of the gap equation may not converge for 
a potential with a strongly repulsive core, as is the case 
with the MF-III. We employ a modified iterative proce- 
dure, which we now describe. 



The starting point is the gap equation with an ultra- 
violet momentum cutoff A ^ Ap, where Ap is of the 
order of the natural (soft) cutoff of the potential. Suc- 
cessive iterations, which generate approximant to the 
gap function from approximant A^*"^-* (? = 1, 2, . . .), are 
determined by 



(27r)3 
A(*-i)(p',A) 



V-''^'{p,p') 

[/(c.+ )-/(c._)]. (17) 



^E} + D(^-^)(p',A) 



The process is initialized by first solving Eq. (|13|l for 
D{pp), where pp is the Fermi momentum, assuming 
the gap function to be a constant. This sets the scale 
of the gap function. The initial approximant for the 
momentum-dependent gap function is then taken as 

^'^^=^){p) = V{pF,p)D{pF). 

There are two iteration loops. The internal loop op- 
erates at fixed A and solves the gap equation H17|l iter- 
atively for i = 1,2,.... The external loop increments 
the cutoff A until the gap becomes insensitive to A, i.e., 
dA{p, A)/dA = 0. The finite range of the potential guar- 
antees that the external loop converges once the entire 
momentum range spanned by the potential is covered. 
Thus, choosing the starting A small enough that the 
strong repulsive core of the potential is eliminated, we 
execute the internal loop by inserting A^'~^^(p) in the 
r.h. side of Eq. H17I) to obtain a new Ai{p) on the l.h. side, 
which in turn is re-inserted in the r.h. side. This proce- 
dure converges rapidly to a momentum-dependent solu- 
tion for the gap equation for A{p, Aj)9{Aj —p), where 6 is 
the step function and the integer j counts the iterations 
in the external loop. 



For the next iteration, the cutoff is incremented to Aj — Aj-i + SA, where SA ^ Aj, and the internal loop is 
iterated until convergence is reached. The two- loop procedure is continued until Aj > Ap, after which the iteration 
is stopped, a final result for A(p) independent of the cutoff having been achieved. Once this process is complete, 
the chemical potential must be updated through Eq. Accordingly, a third loop of iterations seeks convergence 

between the the output gap function and the chemical potential, such that the starting density is reproduced. (This 
loop is not mandatory, since one may choose to work at fixed chemical potential rather that at fixed density.) 

The top panel of Fig. 1 shows the dependence of the gap function on temperature for several densities n, given 
in terms of the ratio / — no/n, where uq = 0.16 fm^"^ is the saturation density of symmetrical nuclear mat- 
ter. The bottom panel shows the associated chemical potentials /i computed self-consistently from Eq. I|14|) . The 
low- and high-temperature asymptotics of the gap function are well described by the BCS relations A(r — > 0) = 
A(0) - [27raA(0)r]i/2exp(-A(0)/r) and A(T ^ T^a) = 3.06^[rc2(Tc2 - T)]!/^, respectively, where Te2 is the critical 
temperature of the phase transition. However, the BCS weak-coupling values a = 1 = P must be replaced with 
a ~ 0.2 and /? ~ 0.9. As a consequence, the ratio of the gap at zero temperature to the critical temperature deviates 
from the familiar BCS result A(0)/Tc2 = 1.76. The deviations from BCS theory are understandable in that (i) the 
system is in the strong-coupling regime, and (ii) the pairing is in a spin-triplet state rather than spin-singlet [i.e., 
D{p)^Aip)]. 

One measure of coupling strength is the ratio A(0)/|/i| of the zero-temperature energy gap to the magnitude of 
the chemical potential. It is seen from Fig. 1 that the strong-coupling regime is realized for / > 40 (i.e. A ^ fj,). At 
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FIG. 1: Dependence of the pairing gap (upper panel) and the chemical potential (lower panel) on temperature for fixed values 
of the ratio / — no/n, where n is the baryon density and no = 0.16 fm~^ is the saturation density of symmetrical nuclear 
matter. Values of the diluteness parameter na^ assume a scattering length a — 5.4 fm. 



f — 20 the system is in a transitional regime (A ~ /i). Another measure of coupling is the diluteness parameter n|ap, 
where a is the scattering length. In agreement with the above criterion, the matter is in the dilute (or strong-coupling) 
regime for / > 40, since this corresponds to na^ = 0.63 < 1 when a is taken as the triplet neutron-proton scattering 
length, 5.4 fm. A signature of the crossover from weak to strong coupling is the change of the sign of the chemical 
potential, which occurs for / w 80 (Fig. slightly below the crossover density between weak- and strong-coupling 
regimes. 

In the limit of vanishing density, / — > oo, the value of the chemical potential at T = tends to fJ,(oo) = —1.1 MeV, 
which is half the binding energy of the deuteron in free space 14]. Indeed, in this limit the gap equation reduces to 
the Schrodinger equation for the two-body bound state, with the chemical potential assuming the role of the energy 
eigenvalue Thus, the BCS condensate of Cooper pairs in the 351 state evolves into a Bose-Einstein condensate 
of deuterons as the system crosses over from the weak- to the strong-coupling regime. The crossover is smooth, taking 
place without change of symmetry of the many-body wave f unction. 



D. Two-body bound states 

Next let us consider temperatures above the critical 
temperature of pair condensation. The two-body T- 
matrix that sums up the particle-particle ladders for a 
system interacting with the potential V obeys the oper- 
ator equation 

T + VGoT + TGoV. (18) 



Since the potential is time-local, the T-matrix depends 
on two time arguments (instead of four in general) , and 
its transformation properties are identical to those of 
the two-point correlation functions discussed above. The 
four-momenta P = {E,P) and p = (e,p) in the center- 
of-mass system are related to their counterparts fci.2 — 
{loi^2, fci.2) in the laboratory system through P = ki + k2 
and p = (ki — k2)/2. 

In the momentum representation, Eq. H18|l takes the 
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form 



T^ip,p';P,E) = V{p,p') + J ^V{p,p") 

xG^{p",P,E)T''{p",p';P,E) (19) 

for the retarded component of the T-matrix. The rele- 
vant two-body Green's function is 



G^(fcl,fc2,S) = 



G>(fci)G>(fc2)-G<(fci)G<(fc2) 

2 E - LUi - UJ2 + IT] 

Q2{ki,k2) 



E - e{ki) - eik2) + nf 



(20) 



where the second relation follows in the quasiparticle ap- 
proximation and we introduce the abbreviation = 
J duj /{2Tr). The two-body phase-space occupation factor 
Q2{ki,k2) = 1 — /(fci) — f{k2), operating in intermedi- 
ate states, allows for propagation of particles and holes, 
thereby incorporating time-reversal invariance. The two- 
body T'^^-matrix has a pole at the energy corresponding 
to the two-body bound state. If Q2 = 1, the pole is ex- 
actly at the binding energy of the deuteron; otherwise the 
pole on the real energy axis determines the binding en- 
ergy of a dimer in the background medium of finite den- 
sity and temperature. Numerical solutions of Eq. H19|) 
will be considered after treating the problem of trimer 
binding in the next section. 



III. THREE-BODY BOUND STATES 

For completeness, we first recapitulate the three-body 
equations 22] that are used in Subsec FlII Bl at finite den- 
sity and temperature to obtain an integral equation for 
the wave function of the three-body bound states. Read- 
ers more concerned with the numerical results can go 
immediately to Subsec. 1111 ij| 



A. Formalism 

A fermionic system supports three-body bound states 
when there exists a non-trivial negative energy solution 
to the homogeneous counterpart of the three-body equa- 
tion 



T = v + vgv = v + vgoT, 



(21) 



for the three-body T-matrix, where V is the three-body 
interaction and Q and Qq are the full and free three-body 
Green's functions. For compactness Eq. (|21|l is written 
in the operator form. If we assume that the particles of 
the system interact via two-body forces, the three-body 
interaction between any three particles (123) reduces to 
a sum of pairwise interactions V = V12 -I- V23 + V13 , where 
Vij is the interaction potential between particles i and j. 
The kernel of Eq. H21|l is not square integrable, since the 



pair potentials introduce momentum-conserving delta- 
functions for the spectator non-interacting particle. As a 
consequence, the iteration series contain singular terms 
(e.g. of type VijGoVij to the lowest order in the interac- 
tion). However, a complete resummation of the ladder 
series in any particular two-body channel ij takes care of 
this problem. 

We decompose the three-body scattering matrix as 
r = r(i) +r(2) +r(3), where 



(22) 



and ijk ~ 123, 231, 312. The channel ij transition oper- 
ators Tij resum the successive iterations with the driving 
term Vij , according to 



(23) 



Eqs. (|22|l and H23|l are combined to eliminate the two- 
body interactions in favor of the channel matrices and 
arrive at a system[3 0| 



(24) 



of three coupled, nonsingular integral equations of Fred- 
holm type II. The Tij-matrices are essentially the two- 
body scattering amplitudes, embedded in the Hilbert 
space of three-body states. Analysis of the bound-state 
problem is based on the homogeneous version of Eq. H22|l . 

The medium modifications encoded in the three-body 
propagator Go become apparent when it is written in the 
momentum representation p^ . 

go{ki,k2,k3,n) = f \G>{h)G>{k2)G>{k3) 

J UJi .UJ2 -,^^3 

i<{u. \n<(h^\n<( 



-G<{ki)G<{k2)G<{k^) 

Q3(fcl,fc2,fc3) 

f2 - e(fci) - e(fc2) - e(fe3) + W 



(25) 



where the ki = {uji,ki) are the particle four- momenta 
and 

Q3(fel,fc2,fc3) - [l-/(fci)][l-/(fc2)][l-/(fc3)] 

- /(fel)/(fc2)/(fc3) (26) 

is the intermediate-state phase-space occupation fac- 
tor for three-particle propagation. In the second line 
of Eq. H25|l . the intermediate-state propagation is con- 
strained to the mass shell within the quasiparticle ap- 
proximation. The momentum space for the three-body 
problem is spanned by the Jacobi four-momenta K — 
ki + kj + kk, Pij = (fci-fcj)/2, an qk = {ki + kj)/3-2kk/3. 
The expressions for particle momenta in terms of Ja- 
cobi coordinates, namely ki — K/3 -I- P13 -I- 92/2, fc2 — 
K/3—prs+q2/2, and k^ = K/3 — q2, are to be substituted 
into Eq. if^ . 
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To obtain the bound-state energy, it is convenient to 
work with the wave function components rather than 
the T^'^ matrices. The total wave function of the three- 
body state is given by the sum 



(27) 



of its three components, which, in analogy to Eq. H22|l . 
obey the homogeneous equations 



(28) 



For identical particles these equations reduce to a single 
equation 

^^'^ = goT23{Pl2P23 + ^13^23)^-^'^ = eoT23PV^'\ (29) 

where P — P12P23 + P13P23 and permutes the indices 
i and j. The total wave function is obtained as ^' = 

(i-hP)V'^^^- 



B. Solving for bound states 

Eq. H29(l can be reduced further to an integral equa- 
tion in two continuous variables by working with states 



diagonal in the angular-momentum basis. 



1, 



\pqa), = \pq{l\)LM{s-)SMs),, 



(30) 



where p and q are the magnitudes of the relative mo- 
menta of the pair {fcj} (indicated by the complementary 
index i), I and A are their associated relative angular- 
momentum quantum numbers, s is their total spin, and 
LMSMg are the orbital and spin quantum numbers of 
the three-body system. The channel T-matrix of the pair 
{fcj} in this basis takes the form 



{pqa\%j\p' q' a') = 5a,a' 



5{q 



-T(p,p\E-qy2m), 



(31) 



while, with to = E — q^ /2m, the two-body T-matrix in 
the two-particle space has the structure 

{plmisms\T{uj)\p'l'm'is'm'^) = 5w5 {p\T{^)\p') 



(32) 



The 



states 



are 



such 



normalized 

iplmismslp'l'm'is'm'^) = p^'^d{p - p')dwSm^„i'Sss'Sr, 



that 



The two-body T-matrix is governed by Eq. (fH^ . rewritten as 



{p\Tico)\p') = {p\V\p') 



dp"p"^ 
47r2 



{p'\V\p") 



Q2{v, q) 



'^-e+(g,p)-e-('7,p)+i?? 



{p"\T{u;)\p') 



(33) 



where Q2{q,p) = (1 ^ /('7/2 + p) — /('7/2 + p)) and e±{q,p) — (e(g/2 ±p)) are averaged over the angle between the 
vectors q and p. The three-body propagator in the \pqa)i basis has the form 



{pqa\go{i^)\p'q'a') = Saa' 



5{p-p')5{q-q') 



Q3iq,p) 



pi 



(34) 



with Q3{q,p) given by the angle average of Eq. H2t)|l . Here the propagator is assumed to be independent of the 
momentum of the three-body system with respect to the background [K ~ 0). Finally, the required expression for 
the permutation operator P in the chosen basis is 



{pqa\P\p'q'a') = 



^ (5(7ri -p) 6(7:2 -p') TT t ' \ 

1 



P' 



[p'f 



(35) 



In this expression, ttI — q'"^ + q^ / A + qq'x and irl 



q'2y4 _|_ qq'x, where x is the angle formed by q and q' , while 



Haa{q,q' ,X) ^^Pn{x) ^ ^ 



ql2+l'2 ^q'y2+l'2 f^nhl[l2l'2 



(36) 



n=0 



h+l2=l I' 



In turn, P„(a;) denotes the Legendre polynomial, and the coefficient s are combinations of 3j and 6j sym- 



bols (31|. The resulting integral equations can be solved by iteration 26] 

Compared to the free-space problem, the three-body equations in the background medium now include two- and 
three-body propagators that account for (i) the suppression of the phase-space available for scattering in intermediate 
two-body states, encoded in the functions Q2 ^ fih) — .f{kj), (ii) the phase-space occupation for the intermediate 
three-body states, encoded in the function Q3 = [I — .f{ki)][l — f{kj)][l — f{kk)] — f{ki)fikj)f{kk), (iii) renormalization 
of the single particle energies e{p) (although the numerical calculations to be reported are carried out with the 
free single-particle spectrum). For small temperatures the quantum degeneracy is large and the first two factors 
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FIG. 2: Dependence of the two-body (Ed) and three-body (Et) binding energies on inverse temperature, for 
fixed values of the ratio / — no/n, where n is the baryon density and no = 0.16 fm~^ is saturation density of 
nuclear matter. For asymptotically large temperature, Ed{oo) = —2.23 MeV and Et{oo) — —7.53 MeV. The 
ratio Et{l3) / Ed{(3) is a universal constant independent of temperature. 



significantly reduce the binding energy of a three-body bound state; at a critical temperature T^a corresponding to 
Et{(3) = 0, the bound state enters the continuum. 



This feature is illustrated in Fig. [3 which shows the tem- 
perature dependence of the two- and three-body bound- 
state energies in dilute nuclear matter for several values 
of the density of the environment. The effect of tempera- 
ture on the two-body bound state is due to the Q2 factor 
only. In analogy to the behavior of the in-medium three- 
body bound state, the binding energy of the two-body 
bound state enters the continuum at a critical tempera- 
ture Tc2, corresponding to the condition Ed{f3) = 0. 

The solutions obtained exhibit a remarkable feature: 
the ratio r/ — Et{l3) / Edifi) is a constant independent of 
temperature. For the chosen potentials, the asymptotic 
free-space values of the binding energies are ^^((O) — 
-7.53 MeV and Ed{Q) = -2.23 MeV; hence 77 = 3.38. 
An alternative definition of the critical temperature for 
trimer extinction is Et{f3!^^) — Ed{(3). This definition 
takes into account the break-up channel t — s- d -|- n of the 
three-body bound state into the two-body bound state 



d and a nuclcon n. The difference between the two def- 
initions is insignificant, because of the property ri{/3) — 
Const. The binding energies of the two and three-body 
bound states, as well as their ratio, depend on the choice 
of the interaction. The universality pointed out above ap- 
plies only to the temperature dependence of these quan- 
tities. It should be tested for other interactions in the 
future. 

Fig. El depicts the normalized three-body bound-state 
wave function for three representative temperatures, as 
a function of the Jacobi momenta p and q. (The tem- 
perature decreases from top to bottom). As the tem- 
perature drops, the wave function becomes increasingly 
localized around the origin in momentum space. Cor- 
respondingly, the radius of the bound state increases in 
r-space, eventually tending to infinity at the transition. 
The wave-function oscillates near the transition temper- 
ature (bottom panel of Fig.|3J). This oscillatory behavior 
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FIG. 3: Wave function of the three-body bound state as a function of the Jacobi momenta p and q defined in 
the text, for / = no/n = 60 and temperatures (top to bottom) T = 60, 10, and 6.6 MeV. 



is a precursor of the transition to the continuum, which 
in the absence of a trimer-trimer interaction is character- 
ized by plane- wave states. 



IV. SUMMARY AND OUTLOOK 

The complexity of the phase diagram of low-density, 
finite-temperature nuclear matter is twofold: (i) The sys- 
tem supports liquid-gas and superfluid phase transitions. 
Thus, depending on the temperature and density, nu- 
clear matter can be in the gaseous, fluid, or superfluid 



state. These features are generic to systems of fermions 
interacting with "van der Waals type" forces that are 
repulsive at short range and attractive at large separa- 
tions, (ii) The attractive component of these forces leads 
to clustering. The clustering pattern is non-universal and 
depends on the form of the attractive force in a particular 
system. 

Fig.Qlcombines results from Sees. Ill and IV in plots of 
the critical temperatures for the superfluid phase transi- 
tion {Tc2) and for extinction of three-body bound-states 
(Tcs), as functions of density. The phase diagram di- 
vides into several distinct regions: (A) The low-density, 
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3 

na 



FIG. 4: Dependence of the critical temperature of extinction 
of three-body bound-states on density (solid line). Trimers 
exist above this critical line. The dependence of the critical 
temperature of the superfluid phase transition is shown by 
the dashed line. The condensate is a weakly-coupled BCS 
superfluid far to the right of the vertical line and is a Bose- 
Einstein condensate of tightly-bound pairs far to the left, with 
a smooth crossover transition in between. 



high-temperature domain is populated by trimers, which 
enter the continuum when the critical line T^sin) is 
crossed from above. (B) The low-temperature and low- 
density domain (na^ <C 1 , / < 40) contains a Bose 
condensate of tightly-bound deuterons. (C) The low- 
temperature, high-density domain features a BCS con- 
densate of weakly-bound Cooper pairs {na^ ^ 1). (D) 
The domain between the two critical lines contains nu- 
cleonic liquid. The phases C and D are characterized 
by broken symmetry associated with the (ipip) conden- 
sate. The transition C^B does not involve symmetry 
changes and is a smooth crossover from the BCS to BE 
condensate. The transitions B-^D and C-^D are second- 
order phase transitions related to the vanishing of the 
condensate along the line Tc2{n). The transition A^D 
can be characterized by an order parameter given by the 
fraction of trimers, which goes to zero at the transition. 
This transition appears to be of second order; however, a 
study of thermodynamics of the transition will be needed 
for confirmation. 

To summarize, we have studied the temperature- 
density phase diagram of dilute isospin-symmetric nu- 
clear matter, which features an isospin-singiet, spin- 
triplet pair condensate at low temperature and a gas 
of trimers (three-body bound-states) at high tempera- 
ture. Solving self-consistently for the gap function and 
the chemical potential, we have quantified the behavior 
of the system in the density-temperature domain where 
the Cooper condensate crosses over to a BEC conden- 



sate of tightly-bound deuterons. A modified iterative 
procedure for solution of the gap equation with a "run- 
ning" cutoff has been designed and implemented for this 
purpose. The method is found to be effective in treat- 
ing the strong repulsive core of the potential. Solution 
of the three-body integral equations for bound states in 
the background medium furnishes us with density and 
temperature-dependent binding energies of trimers in di- 
lute nuclear matter. Numerically, we find that the ratio 
of the temperature-dependent binding energies of dimers 
and trimers is independent of the temperature and den- 
sity and can be determined from its value in free space. 
For large degeneracies, the binding energy of trimers is 
small and vanishes at a critical temperature Tc3 which is 
larger than the critical temperature of superfluid phase 
transition Tc2- These critical lines Tc3(n) and Tc2{n) sep- 
arate the phase diagram into distinct regions. At the 
high-temperature/low-density end the system is popu- 
lated by trimers, whereas in the low-temperature/high- 
density regime the system supports a condensate of 
neutron-proton Cooper pairs, which crosses over to a 
Bose-Einstein condensate of deuterons as the density de- 
creases. 

We have omitted a number of correlation effects in our 
discussion, (i) The momentum and energy dependent 
self-energies renormalize the spectrum of quasiparticles 
in nuclear matter. While these renormalizations are im- 
portant at densities comparable to the nuclear saturation 
density, it is likely that their influence on the physics at 
much smaller densities (/ < 40) is minor, (ii) It is known 
that the screening of interactions, associated with long- 
wave length density and spin fluctuations, are important 
for dilute fermions; e. g. they change the value of the 
pairing gap by a factor of 0.45 jl5, 32, 33, 3^j. A task for 
the future is to incorporate the screening of interactions 
in the bound state calculations discussed above. 

Clearly, the true phase diagram of dilute and finite- 
temperature nuclear matter is richer than that derived 
from the present study. Modifications are likely to origi- 
nate from clusters with A = 4 and their Bose-Einstein 
condensation (35j . Still, we expect that although the 
phase diagram shown in Fig. ^ has been obtained for 
a particular form of the two-body interaction, its quali- 
tative features should be universal for saturating systems 
of attractive fermions that possess two- and three-body 
bound states in free space. 
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